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loss 0.035 < < 0.095 and 4-momentum transfer squared \tp\ < 1.0 GeV^, the proton fractional 
momentum loss to the system X was evaluated from the momenta of the particles comprising X. 
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exchange. The fraction of events with ^p < 0.02 is found to be 0.194 ± 0.001 (stat) ± 0.012 (syst). 

PACS numbers: 11.55.Jy, 12.40.Nn 



The success of perturbative Quantum Chromo- 
Dynamics (QCD) in describing strong interactions at 
high transverse momentum transfers rests on the fac- 
torization theorem, which allows hadronic cross sec- 
tions to be expressed in terms of parton-level cross 
sections convoluted with uniquely defined hadron par- 
ton densities. It is therefore not surprising that the re- 
cently reported breakdown of factorization in diffrac- 
tive dijet production 0, a process containing both a 
hard scattering and the characteristic rapidity gap sig- 
nature of diffraction, has attracted considerable theo- 
retical attention. Rapidity gaps, defined as regions of 
pseudorapidity devoid of particles, are presumed 
to be formed in diffractive events by the exchange of 
Pomerons {IP), which in QCD correspond to entities 
of gluons and/or quarks with the quantum numbers 
of the vacuum Q (see Fig. 1). The breakdown of fac- 
torization in diffraction is generally attributed to a 
suppression of the cross section resulting from addi- 
tional partonic interactions within a diffractive event 
that spoil the rapidity gap signature 0, 0| ■ In pro- 
cesses with two rapidity gaps, as in that with two for- 
ward gaps traditionally referred to as double Pomeron 




FIG. 1: Diagrams and event topologies for (a) single 
diffraction, p -\- p —> p -\- X , and (b) double Pomeron (_P) 
exchange, p-\-p—>p + X + Y; the shaded areas represent 
pseudorapidity regions of particle production. 



exchange (DPE), shown in Fig. lb, one might then 
expect that either both gaps survive or are simultane- 
ously spoiled, leading to a largely non-suppressed ratio 
of two-gap to one-gap rates ,6] . Such a scenario could 
explain our finding that the ratio of the rates of DPE 
to single diffractive (SD) dijet production is about five 
times larger than that of SD to non-diffractive (ND) 
dijet production 0. 



4 



Since rapidity gap formation is a non-perturbative 
phenomenon, soft (low transverse momentum) difFrac- 
tive cross sections would be expected to exhibit a simi- 
lar behavior. Indeed, the SD pp cross section has been 
found to be suppressed at high energies by a factor of 
^ 10 relative to extrapolations from lower ener gy d ata 
based on Regge theory and factorization Hlill^. In 
this Letter, we present a measurement of the ratio 
of the inclusive DPE to SD cross sections in pp colli- 
sions at a/s = 1800 GeV and compare our results with 
previous measurements |11| and with predictions from 
Regge theory and various theoretical models proposed 
to account for the breakdown of Regge factorization 
in SD. Our measurement severely constrains the avail- 
able models, paving the way towards a more compre- 
hensive understanding of the physics of rapidity gaps. 

The components of the Collider Detector at Fermi- 
lab (CDF) relevant to this study are the Roman Pot 
Spectrometer (RPS) 0, used to detect leading an- 
tiprotons, and the calorimeters and beam-beam coun- 
ters (BBC) 0, used to detect the particles from pro- 
ton dissociation. The RPS measures the fractional 
momentum loss and 4-momentum transfer squared 
tp of the antiproton with resolutions 5£p = ±1.0x 10"^ 
and Stp = ±0.07 GeV^, respectively [ij. The calorime- 
ters have projective tower geometry and cover the re- 
gions |?7| < 1.1 (central), 1.1 < I77I < 2.4 (plug), and 
2.2 < \r]\ < 4.2 (forward). The Ary x Acj) tower dimen- 
sions are approximately 0.1 x 15° for the central and 
0.1 X 5° for the plug and forward calorimeters. The 
BBC consist of two arrays of eight vertical and eight 
horizontal scintillation counters perpendicular to the 
beam line at z = ±6 m, BBCp and BBCp, covering 
approximately the region 3.2 < \ri\ < 5.9 in four 77- 
segments of width A77 « 0.7. 

The present study is based on our ^/s =1800 GeV 
inclusive SD data sample 0. The events were col- 
lected in the 1995-96 Tevatron Run IC by trigger- 
ing on an antiproton detected in the RPS. Offline 
cuts were applied requiring a reconstructed track in 
the RPS, no more than one reconstructed vertex in 
the CDF detector within a distance |zvtx| < 60 cm 
from the nominal beam-beam interaction point along 
the beam direction, and a BBCp multiplicity of < 6. 
These cuts remove overlap events due to multiple in- 
teractions in the same beam-beam crossing, compris- 
ing 4% of the inclusive SD data sample as estimated 
by the instantaneous luminosity. 

Experimentally, we study the DPE process p + p ^ 
p + X + Y , where Y is either a leading proton or 
a low-mass proton dissociation system which escapes 
undetected through the beam-pipe on the proton side; 
the mass-squared of the system Y is estimated to be 
My ^ 8 GeV^. The procedure we follow to identify 
and measure the DPE signal in these data is to select 
an event sample with {^p, tp) within a certain region 
and measure the fractional longitudinal momentum 



£,p of the proton transferred to the system X using 
the equation 

1 " 

^ i—l 

where and 77* are the transverse energy and pseu- 
dorapidity of a particle and the sum is carried out 
over all particles excluding the proton or the parti- 
cles associated with the system Y. DPE events are 
expected to appear in the low region, in contrast 
to SD events for which w 1. In practice, not all 
particles of the system X are included in evaluating 
Eq. (1) because (a) CDF does not provide full coverage 
and (b) particles depositing energy in the calorimeters 
below the energy thresholds used to reject noise are 
excluded. This issue is addressed by applying appro- 
priate correction factors and by calibrating formula 
(1) on the antiproton side by directly comparing the 
value of obtained by this method with that mea- 
sured by the RPS, ^p^^, as discussed below. 

To evaluate we use calorimeter towers and BBC 
hits. The tower energy thresholds used, chosen to 
lie comfortably above noise level, are Et — 0.3 GeV 
for the central, Et = 0.2 GeV for the plug, and 
E — 1.5 GeV for the forward calorimeters; at the 
calorimeter interface near |?7| ~ 2.4 a threshold of 
Et = 0.275 GeV was used. These values are based 
on test-beam calibrations of the calorimeters 12] and 
must be multiplied by an Ty-dependent factor /^j^ (of 
average value (/b^) = 1.6) to obtain the true Et at 
low energies [l3| • To account for particles below tower 
threshold, the calorimeter contribution to is mul- 
tiplied by /thr = 1.54. This factor is obtained from a 
Monte Carlo (MC) simulation in which the same tower 
thresholds are used as in the data after dividing the 
generated particle energy by Jet- The Monte Carlo 
simulation is based on the single diffractive genera- 
tor described in |3| and references therein, adapted to 
double Pomeron exchange. For each BBC hit we use 
77 and Et values randomly chosen from a flat rj dis- 
tribution over the hit BBC ?7-segment and from the 
shape of the Et distribution expected from the MC 
simulation, respectively. The BBC contribution to 
is then weighted by a factor of 3/2 to account for neu- 
tral particles, which are undetected by the BBC, and 
by an additional factor of 3/4 to account for the over- 
lap regions among the four scintillation counters of 
each BBC segment. Hits in the outer ?7-segments, 
3.2 < |ry| < 3.9, which overlap with the forward 
calorimeters, are ignored. The BBC contribution to 
is less than 10% in the region of -4 < logf^ < -2 
and increases to 60% at log^^ = — 5 and log^^ = — 1. 

The method of measuring ^ using Eq. (1) is cal- 
ibrated on the antiproton side by evaluating = 

J27=i ETe~^ (excluding the antiproton from the 
sum) and comparing its value with that measured by 
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the RPS. The data are divided into bins of A^^^^'^' = 
0.01, and the £,p values obtained for each bin are fit- 
ted with a Landau distribution. Figure 2a shows, as 
an example, the data and fit for 0.05 < £,p^^ < 0.06. 
The ratio of width to peak position is « 0.6 over the 
entire £_p region of our data sample. The enhance- 
ment in the small £,p region is caused by a downward 
shift in in low multiplicity events due to "loss" of 
particles with energy under tower threshold. Within 
the region 0.01 < £,p^^ < 0.1, an approximately lin- 
ear relationship is observed between the median value 
of and ^^P^. A fit with « CC|^^^ yields 
C — 0.95, in close agreement with the expected value 
C — 1. A fit in which C = 1 and (/b^) is varied with 
/coir = {Jet ) ^ /thr treated as a free parameter yields 
/coir — 2.7. In Fig. 2b an error of ±5% is used in 
all data points to yield x^/d.o.f. = 1 for this fit. In 
extracting results, we use /corr = 2.7 and assign a con- 
servative ±10% error to C (twice the error obtained 
from the fit) to account for other possible systematic 
uncertainties. 
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FIG. 2: (a) Distribution of antiproton fractional mo- 
mentum loss measured from calorimeter and beam- 
beam counter information for events in wliicli the (,p^^ 
value measured by the Roman Pot Spectrometer is within 
0.05 < ^^^'^ < 0.06; the sohd line is a Landau fit. (b) 



Median values ^p obtained from Landau fits to data in 
different bi 
ship is observed. 



different £,p^^ bins plotted versus ^p^^\ a linear relation- 



The DPE signal is evaluated for events with an- 
tiproton S^p and tp within 0.035 < < 0.095 and 
\tp\ < 1.0 GeV^, where the RPS acceptance is larger 
than « 30% 0. The total number of inclusive SD 
events in this region is 568K. The caUbrated dis- 
tribution is compared in Fig. 3 with a two-component 
MC simulation that includes SD and DPE. The shape 
of the input ^p distribution in the MC simulation for 
DPE is based on a triple-Pomeron term on the pro- 
ton side using a Pomeron intercept a]p{0) = 1 + e 
with e — 0.104, as determined from a global fit to 



p{p) total cross section data 0|. DPE events were 
generated for < 0.1. The DPE and SD MC gener- 
ated events are independently normalized to the data 
points in the regions 4 x 10~^ < ^p < 10~^ and 
0.02 < < 1, respectively. The SD events appear 
as a broad peak around = 1, which falls exponen- 
tially as ^p decreases. The DPE events appear as a 
flattening of the distribution on the low side and 
represent the dominant contribution for < 0.02. 
The wavy shape of the data distribution in the DPE 
region is due to the ry-dependent calorimeter tower 
energy thresholds used and is reproduced by the MC 
simulation. At low both data and MC simula- 
tion extend down to and below the kinematic limit of 
Cp,min = M^/(sCp,min) ~ 10~^, whcre Mq is the low- 
est mass for DPE excitation after threshold turn-on 
effects set in, taken to be 1 GeV. The events below 
the kinematic limit are due to the downward fluc- 
tuations of in low multiplicity events mentioned 
above. The agreement between data and MC simu- 
lation in the region of < 0.02 shows that Regge 
factorization is successful in describing the shape of 
the ^ distribution in DPE using the Pomeron inter- 
cept determined in ^14]. The ratio of the number of 
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FIG. 3: Distribution of proton fractional momentum loss 
, measured from calorimeter and beam-beam counter 
information, for events with a leading antiproton of 
0.035 < < 0.095 and \tp\ < 1.0 GeV^ the curves 

are from a Monte Carlo simulation of SD (dotted), DPE 
(dashed) and total (solid) contributions normalized to the 
data points; the DPE events were generated for < 0.1. 



events within < 0.02 to the total number of events 
is 0.202 ± 0.001 (stat). After correcting for smearing 
effects caused by the resolution, the ratio becomes 
rDpe = 0.194 ± 0.001 (stat) ± 0.012 (syst), where 
the systematic error is from the uncertainties due to 
calibration (±0.003), smearing (±0.008), and 
low enhancement (±0.008, see Fig. 2a) added in 
quadrature. 
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Neglecting Reggeon contributions, the DPE/SD ra- 
tio is given in Regge theory by ^\ 



R 



DPE\ 
SD \ip 



0.02 



K/3'^{tp)dtpd^p 



(2) 



where k is the ratio of the triple-Pomeron cou- 
pHng g(tp) to the Pomeron-proton couphng (3{tp) and 
a{t) = a{0) + a't is the Pomeron trajectory. Using 
K = 0.170 ±0.017, P{tp) = /3(0)e4-6*P, p^ioyiGn = 
0.86 GeV-2 and a{t) = 1.104 + 0.25t jl3 yields 
R§^^ (Regge) = 0.36 ±0.04. This prediction is larger 
than the measured values by a factor of 1.9±0.2. How- 
ever, this discrepancy from the factorization expecta- 
tion of unity is small compared to the 0(10) discrep- 
ancy observed in SD 9J. Thus, this result confirms the 
conjecture made in the introduction that the forma- 
tion of a rapidity gap within the rapidity space cov- 
ered by the diffraction dissociation products in events 
with a leading (anti)proton would be largely non- 
suppressed. A similar conclusion has been reached 
by the UA8 Collaboration from a study of DPE pro- 
duction in pp collisions at ^/s =630 GeV at the CERN 
SppS collider [llj. Changes in the predicted two-gap 
to one-gap ratio due to contamination of the DPE sig- 
nal with proton fragmentation events are estimated to 
be ^ 15% and therefore are not expected to alter this 
conclusion. 

Phenomenological models proposed to account for 
the breakdown of Regge factorization in SD may be 
divided into two broadly defined classes: (a) those 
attributing the violation either to "damping" of the 
cross section at small ^ or to a decrease of the 
Pomeron intercept at low ^ jl6l | or at high ener- 
gies and (b) those in which the overall normal- 
ization decreases with increasing energy but the shape 
of the ^ distribution remains practically 4, 5, 1^ or 
entirely 0, 0| unchanged. The models of class (a) 
predict a ^p distribution different from that expected 
from SD and are disfavored by the shape of the distri- 
bution presented in Fig. 3, which behaves as l/^"("' 
down to the kinematic limit of ^p,min ~ 10~^. Of 
the class (b) models, three have reported predictions 
for both SD and DPE: the eikonal model the 
Pomeron flux renormalization modellsl , and the gap 
probability renormalization model l6f|. The eikonal 
model, in which "screening corrections" to the Regge 
amplitude are calculated using an eikonal approach, 
yields suppression factors of 0.369 and 0.309 for SD 
and DPE, respectively; although the DPE/SD ratio 
is relatively non-suppressed, in close agreement with 
our result, the suppression for SD is underestimated 
by a factor of ~ 3. The Pomeron flux renormaliza- 
tion model, in which the Regge theory Pomeron flux 
factor is renormalized to unity for Pomerons emit- 
ted by the p in SD or DPE and independently by 
the p in DPE, yields the correct suppression factor 
for SD, but predicts a DPE/SD ratio smaller than 



the measured value by a factor of 4.7 ±0.6 [9j. Fi- 
nally, in the gap probability renormalization model, 
in which the SD and DPE cross sections are expressed 
in terms of the variables M|- and Ar/ — Arjp + Ar/p, 
where Arji = — Inf^, the predicted DPE/SD ratio is 
0.21 ± 0.02 1^, in good agreement with our measured 
value of 0.194 ± 0.001 ± 0.012. These predictions do 
not include possible effects from Reggeon exchange or 
contributions from proton fragmentation. 

In summary, we have studied the double Pomeron 
exchange (DPE) process p + p ^ p + X + Y, where 
y is a leading proton or a proton-dissociation system 
of mass-squared My ^ 8 GeV^, by measuring the 
fractional longitudinal momentum loss of the proton 
to the system X, £^p , in events with a leading an- 
tiproton of 0.035 < ^,5 < 0.095 and \tp\ < 1.0 GeV^ 
produced in pp collisions at ^/s =1800 GeV. Events in 
the region < 0.02 follow a distribution of the form 

^ ^/^p^^\ where a(0) is the Pomeron intercept, and 
are attributed to DPE production. The ratio of the 
number of DPE events in this region to the total num- 
ber of SD events is found to be 0.194 ± 0.001 ± 0.012. 
This value is lower than the prediction based on Regge 
factorization by a factor of 1.9 ± 0.2, which is rel- 
atively small compared to the suppression factor of 
0(10) observed in SD indicating that the forma- 
tion of a second rapidity gap in a SD event is rela- 
tively non-suppressed. Among models proposed to ex- 
plain the suppression of the SD cross section at high 
energies, our results favor those in which the Regge 
based shapes of the SD and DPE distributions remain 
unchanged and only the overall normalization is sup- 
pressedH Ilia, 1,18). 
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